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Gamma-ray bursts (GRBs) are short and intense emission of soft 7— rays, which have fascinated 
astronomers and astrophysicists since their unexpected discovery in 1960s. The X-ray /optical/radio 
afterglow observations confirm the cosmological origin of GRBs, support the fireball model, and 
imply a long-activity of the central engine. The high energy 7— ray emission (> 20 MeV) from 
GRBs is particularly important because they shed some lights on the radiation mechanisms and can 
help us to constrain the physical processes giving rise to the early afterglows. In this work, we review 
observational and theoretical studies of the high energy emission from GRBs. Special attention is 
given to the expected high energy emission signatures accompanying the canonical early-time X-ray 
afterglow that was observed by the Swift X-ray Telescope. We also discuss the detection prospect 
of the upcoming GLAST satellite and the current ground-based Cerenkov detectors. 

PACS numbers: 95.30.Gv, 95.85.Pw, 98.70.Rz 



I. INTRODUCTION 

Gamma-ray bursts (GRBs) are brief intense flashes of 
soft (0.01 — 1 MeV) 7-rays that are detected once or twice 
a day for a BATSE-like detector. GRBs were serendip- 
itously discovered by Vela satellites in late 1960s, and 
were first pubHcly reported by Klebesadel et al. in 1973 
The observed bursts arrive from apparently random 
directions in the sky and they last between tens of mil- 
lisecond and thousands of seconds [2]. Their physical 
origin has been debated for a long time mainly due to 
the lack of an exact position and a reliable estimate of 
the distance to us. In 1997, several GRBs were rapidly 
and accurately localized by the Italian-Dutch BeppoSAX 
satellite, leading to the discovery of their X-ray, optical 
and radio counterparts, and their redshifts [1, 0, B Q- 
The cosmological origin of most, if not all GRBs, was con- 
firmed. The leading interpretation of the data is the cos- 
mological fireball model 

0, i, i, [M [H E E Q d, [3 , 

in which the prompt 7— rays are powered by the baryon- 
rich (or Poynting flux) relativistic jets ejected from the 
central engine with variable Lorentz factors (i.e., internal 
shocks) while the afterglows are produced by the interac- 
tion between the outflow material and the medium (i.e., 
external shocks). 

The most widely discussed radiation mechanisms in- 
clude synchrotron emission and inverse Compton (IC) 
scattering. Both can produce electromagnetic emission in 
a very wide energy range, i.e., radio to hard 7— rays (GeV 
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or even TeV) . Like in Active galactic nucleus ( AGNs) we 
expect that the IC process will give rise to a high energy 
component that will be emitted along with the prompt 
sub-MeV photons and the afterglow radio/optical/X- 
ray emission, as detected in dozens of GRBs by the 
Compton Gamma Ray Observatory (CGRO) satellite in 
1991-2000 [13, [H, [li, H [m, [12 [2|. For example if 
the typical synchrotron frequency of the prompt emis- 
sion is 100 keV and the electrons have a Lorentz factor 
of 500 (in internal shocks) we expect IC emission peaks 
at ~ 20 GeV. Similarly, inverse Compton of reverse shock 
photons with ~ 1 eV by forwards shock electrons with a 
Lorentz factor of 10^ will result in an IC component of 
- 100 MeV. 

The Large Area Telescope (LAT) onboard the 
Gamma Ray Large Area S pace Telescope (GLAST; see 
'http : / /glast . gsfc . nasa. gov /) , to be launched soon, is ex- 
pected to enhance high energy detection rate significantly 
because of its larger effective area than that of EGRET. 
For a bright burst at a redshift z ~ 1, LAT may col- 
lect 10 tens-MeV to GeV afterglow photons The 
estimate of the detectability of the prompt high energy 
7— rays is more difficult because the physical parame- 
ters involved in the internal shocks are still poorly con- 
strained. Regardless of the uncertainties, preliminary 
calculations suggest a promising detection prospect for 
LAT Hi, [2^. 

The high energy emission from GRBs can help us to 
better understand the physical composition of the out- 
flow, the radiation mechanisms, and the underlying phys- 
ical processes shaping the early afterglow. Such goals, of 
course, are hard to achieve because of the rarity of the 
high energy photons. However, with LAT, ^ 10'^ high en- 
ergy photons could be detected from an extremely bright 
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burst (for example, GRB 940217, GRB 030329 and GRB 
080319B) and these can be used to constrain the models. 
With such a hope, we present in this work an overview 
of the theoretical studies of high energy emission from 
GRBs. 

The structure of this review is as follows. We first dis- 
cuss the observational aspects of high energy emission of 
GRBs and afterglows in Section |TT1 and then the physi- 
cal processes in Section Hill We discuss the high energy 
emission processes in GRBs and afterglows, the interpre- 
tations of available high energy observations and possible 
progresses in the next decade in Sections IIVIIVII respec- 
tively. 



II. OBSERVATIONS 

We begin with a short review of the observations of 
high energy emission from GRBs and afterglows. We 
divide this section into four parts, beginning with an in- 
troduction of the detectors. We then have a short discus- 
sion of the cosmic absorption of high energy 7— rays, that 
plays a crucial role in the detection prospects above 50 
GeV. We continue with the prompt high energy emission- 
the GRB itself but in energies above 20MeV and proper- 
ties of high energy afterglows. 



A. Detectors 

Space telescopes. Among space telescopes 
dedicated to high energy 7— ray astrophysics, 
three are particularly interesting for GRB peo- 
ple, including the Energetic Gamma Ray Ex- 
periment Telescope (EGRET) onboard CGRO 
(see http : / /heasarc . gsfc . nas a. gov /docs /cgro /egret / ) , 
Gamma-Ray Imaging Detector (GRID) onboard AGILE 
(see http:/ /agile.rm.iasf.cnr.it/), and the upcoming LAT 
onboard GLAST satellite. EGRET and GRID have a 
similar peak effective area ~ 1000 cm^. That of LAT, 
however, is much larger. We compare them in Tabic [H 
With a higher sensitivity, LAT is expect to detect high 
energy photons from GRBs much more frequently than 
both EGRET and GRID. The Burst Monitor (GBM) 
onboard GLAST is sensitive to X-rays and gamma 
rays with energies between 8 keV and 25 MeV. The 
combination of GBM and LAT provides a powerful 
tool for studying gamma-ray bursts, particularly for 
time-resolved spectral studies over a very broad energy 
band. 

The effective area of the LAT as a function of photon 
energy is shown in the upper panel of FiglT] For photons 
below 100 MeV, the effective area of LAT is small, which 
limits the detection prospect of the MeV photons. GBM 
won't help in this aspect because of its rather small area 
~ 126 cm^. The high energy photons are much less nu- 
merous than the keV-MeV photons because of the large 
energy each photon carries. So even in the most opti- 
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FIG. 1: Upper panel: On-axis effective areas of LAT (from 

http: / /www- glast .slac.stanford.edu / software /IS / glast _ lat _ performance.h 
Lower panel: Effective area of MAGIC for three different 
zenith angles (from [29" 



mistic extreme cases, the number of high energy photons 
detected by GLAST-like satellites is not expected to be 
much more than ^10''. 

Ground-based telescopes. Ground-based high 
energy detectors have very large effective areas 
10* — 10^ but work in the energy range 

of tens GeV to 100 TeV. There are two kinds 
of Cherenkov telescope s: water Cherenkov tele- 
scopes like Milagro fcttp: //www.lanl.gov/milagro/ 1 
and atmospheric Cherenkov telescopes, such 
as MAGIC, H.E.S.S., Whipple, Cangaroo-III 
( http : / / wwwmagic . mppmu.mpg.de /links / ) and 
VERITAS (http://veritas.sao.arizona.edu/ ), and 
their advanced generations, like MAGIC-II 
(http: / / wwwmagic. mpp mu.mpg.de/introduction/magic2[ 
html) and H.E.S.S.-II. 

Milagro is a TeV gamma-ray detector locating in 
northern New Mexico operating in the energy band > 100 
GeV. It uses the water Cherenkov technique to detect ex- 
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TABLE I: LAT Specifications and Performance Compared witli EGRET (from|http:/ /glast. gsfc.nasa.gov/science/instruments/tablel-l.html 
and GRID (from 

Quantity LAT EGRET GRID 

Performance period (year) 2008-? 1991-2000 2006-? 

Energy Range 0.02- 300 GeV 0.02 - 30 GeV 0.03-50GeV 
Peak Effective Area ~ lO'' cm^ 1500 cm^ 700 cm^ 

Field of View 2.5 sr 0.5 sr 3 sr 

Angular Resolution < 3.5° at O.lGeV 5.8° at O.lGeV 4.7° at O.lGeV 



tensive air-showers produced by very high energy (VHE) 
gamma rays as they interact with the Earth's atmo- 
sphere. The Milagro field of view is ~ 2 sr and duty 
cycle is > 90%. The effective area is a function of zenith 
angle and ranges from 50m^ at 100 GeV to 10^ m^ at 10 
TeV [23. 

Among the atmospheric Cerenkov telescopes, MAGIC 
may be best suited for the detection of the prompt emis- 
sion of GRBs, because of its low energy threshold, its 
large effective area, and in particular, its capability for 
fast slewing \2§\. The low trigger threshold, currently 50 
GeV at zenith, should allow the detection of GRBs even 
at a redshift ^ ~ 1, as lower energy radiation can effec- 
tively reach Earth without interacting much with the dif- 
fusive infrared background (see section Hi Bl for details). 
Moreover, in its fast slewing mode, MAGIG can be repo- 
sitioned within 30 s to any position on the sky; in case 
of a target-of-opportunity alert by GGN, an automated 
procedure takes only few seconds to terminate any pend- 
ing observation, validate the incoming signal, and start 
slewing toward the GRB position. So far, the maximal 
repositioning time has been ^ 100 s j2^. In its current 
configuration, the MAGIC photomultiplier camera has a 
field of view of 2.0° diameter and a peak collection area 
for 7— rays of the order of 10^ m^ (see the lower panel of 
Fig[T]) [291]. 

The H.E.S.S. array is a system of four 13m-diameter 
lACTs located at 1 800 m above sea level. The effective 
collection area increases from ~ lO'^m^ at 100 GeV to 
more than lO^m^ at 1 TeV, for observations at zenith 
angles of < 20 deg |30|] . The slew rate of the array is ~ 
lOO'^ per minute, enabling it to point to any sky position 
within 2 minutes j30l |. 

VERITAS is a new major ground-based gamma-ray 
observatory with an array of four 12m optical reflectors 
for gamma-ray astronomy in the energy range of 50 GeV 
— 50 TeV band (with maximum sensitivity from 100 GeV 
to 10 TeV). The telescope design is based on the design 
of the existing 10m gamma-ray telescope of the Whipple 
Observatory. In the energy range from 100 GeV to 30 
TeV, VERITAS'S effective area rises from around 3 x 
lO'^ m^ to well over 10^ m^ and its energy resolution is 
10-20% [3l'|. Its slew rate is similar to that of H.E.S.S.. 
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FIG. 2: Optical depth of the universe to 7-rays from interac- 
tions with photons of the intergalactic background light (IBL) 
and CMB for 7-rays having energies up to 100 TeV. The solid 
lines are for the fast evolution IBL cases, and the dashed lines 
are for the baseline IBL cases. From bottom to top the curves 
correspond to redshift 2 = (0.03, 0.117, 0.2, 0.5, 1, 2, 3, 5), 
respectively (from [34] ) . 

B. Cosmic Attenuation 

As 7— rays with an energy above tens GeV travel to- 
ward the observer, they are absorbed due to the interac- 
tions with the diffuse infrared background [12, HI] . The 
probability of a high energy 7— ray with an energy to 
reach earth without being absorbed is exp(— t), where t 
can be parameterized by 

logT = Ax"^ + Bx^ + Cx^ + Dx + E, 

X = log(e-y/leV) and the coefficients, that approximate 
the lines in Figl^l have been presented in Table HIl 'si'. '3?] . 
For z ^ \ and ~ 100 GeV, we have r ^ 6, which 
limits the detection prospect of the VHE emission from 
GRBs. The spectral energy distribution of intergalactic 
background light (IBL), particularly at redshift z > 0.2, 
is uncertain. So is the estimate of r. In some models, for 
2: ~ 1 and ~ 100 GeV, a r ~ 1 is predicted [36[ (their 
recent estimate, however, gives a r ~ 3.3). If correct, 
the detection prospect of VHE photons from GRBs will 
be more promising. The VHE detection of GRBs will 
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TABLE II: Coefficients for the Baseline IBL (upper low) and 
Fast Evolution IBL Fits (lower low). The parametric approx- 
imation holds for 10"^ < r < 10^ and < 2 TeV for all 
redshifts but also up to 10 TeV for redshifts z < 1 (from 
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A 


B 


C 


D 


E 


0.03 


-0.020228 


1.28458 


-29.1498 


285.131 


-1024.64 




-0.020753 


1.31035 


-29.6157 


288.807 


-1035.21 


0.117 


0.010677 


-0.238895 


-1.004 


54.1465 


-313.486 




0.022352 


-0.796354 


8.95845 


-24.8304 


-79.0409 


0.2 


0.0251369 


-0.932664 


11.4876 


-45.9286 


-12.1116 




0.0258699 


-0.960562 


11.8614 


-47.9214 


-8.90869 


0.5 


-0.0221285 


1.31079 


-28.2156 


264.368 


-914.546 




0.0241367 


-0.912879 


11.7893 


-54.9018 


39.2521 


1.0 


-0.175348 


8.42014 


-151.421 


1209.13 


-3617.51 




-0.210116 


10.0006 


-178.308 


1412.01 


-4190.38 


2.0 


-0.311617 


14.5034 


-252.81 


1956.45 


-5671.36 




-0.397521 


18.3389 


-316.916 


2431.84 


-6991.04 


3.0 


-0.34995 


16.0968 


-277.315 


2121.16 


-6077.41 




-0.344304 


15.8698 


-273.942 


2099.29 


-6025.38 


5.0 


-0.321182 


14.6436 


-250.109 


1897.00 


-5390.55 




-0.28918 


13.2673 


-227.968 


1739.11 


-4969.32 



provide a significant source s amp le for studies of IBL as 
a function of look-back time [33 ■ 



C. Observation of high energy prompt emission 

EGRET had detected more than 30 GRBs with high 

energy photon emission [13, mM US HH, 113,111. Some 

of these > 30 MeV photons are simultaneous with the 
keV— MeV emission, i.e, they are prompt high energy 
emission. 

GRB 940217 is a good example [2l|. This burst 
was detected by the Compton Telescope (COMPTEL), 
EGRET, Burst And Transient Source (BATSE) and 
Ulysses. The blue line in Figl3] is the Ulysses 25—150 
keV light curve. The prompt soft 7— ray emission was 
clearly visible in a timescale of ~ 180 s (i.e, phase-1). 
Simultaneously, 10 photons ranging from 40 MeV to 3.4 
GeV were recorded. The count rate of these photons is 
much higher than that of phase-2, the high energy after- 
glow. 

Milagrito observation of GRB 970417A at energies 
above ~ 0.1 TeV hinted at a distinct higher-energy com- 
ponent (at 3(7 level) , but lacked energy resolution to pro- 
vide a spectrum [39| . The excess had a chance probability 
of 2.8 X 10~^ of being a fluctuation of the background (one 
detection in 54 bursts). Milagro observed more than 50 
GRBs but got null results. The Milagro 99% confidence 
upper limit on the 0.2 — 20 TeV fluence ranges from 10~^ 
to 10"'^ erg cm~^ [2^. Null results at lower energies are 
also reported by MAGIC observations f40| . 

The null results at > 100 GeV are not surprising be- 
cause of the huge absorption of such high energy 7— rays 



by the fireball (see eq. ((50)) below) and by the diffusive in- 
frared background, as already mentioned in section fllB I 



D. Observation of high energy afterglow 

In several GRBs, photons above 30 MeV have arrived 
after the end of the prompt keV— MeV emission. These 
are classified as high energy afterglow. It is very in- 
teresting to note that the afterglow emission has been 
first detected in this energy range rather than in X- 
ray/optical/radio bands even though in low energy bands 
the photons are more abundant by a factor of 10^ — 10^. 

The two best known high energy afterglows are those 
associated with GRB 940217 and GRB 941017. 

GRB 940217: The prompt soft and hard 7— ray emis- 
sion has been described in the last subsection. As the 
25—150 keV emission ceased, the hard 7— ray emission 
did not and lasted longer than 5400 s, including an 18 
GeV photon that arrived about an hour after the trigger 
(see the green circles in phase-2 of Fig|31). In total 18 
high energy photons have been detected. The total num- 
ber could have been higher (probably around 100) if the 
source was not occulted by the earth for ~ 3700 s after 
the burst. Note that the 18 GeV photon was observed 
after the satellite came out from the earth occultation. 
This high energy afterglow is also characterized by: (a) 
The count rate of high energy photons seemed to be con- 
stant; (b) Except of one photon with an extremely high 
energy ^ 18 GeV, the energy of the others is nearly a 
constant, i.e., ^ 100 MeV. 

GRB 941017: GRB941017 is one of highest fluence 
bursts observed by BATSE in its 9-yr lifetime. Ninety 
per cent of the flux observed by BATSE occurred in a 
time interval of 77 s. The high-energy component car- 
ried at least 3 times more energy than the lower energy 
component and it lasted about 3 times longer [H. As 
shown in FigUl there are two additional amazing obser- 
vation facts: (1) While the soft 7— ray emission became 
weaker and weaker and disappeared, both the spectrum 
and the flux of the hard 7— ray emission (up to an energy 
> 200 MeV) were almost constant over a timescale of ^ 
200 s; (2) The spectrum of the high energy emission com- 
ponent is rather hard, F^, oc v'^ where v is the observed 
frequency of the photon. The peak energy of the hard 
component is likely to be above 200 MeV. 

Some ground-based Cherenkov detectors have been 
used to observe the VHE afterglow emission of GRBs. 
After several years' of search, no evidence for afterglow 
photons at such high energies has been found [29,, ^] . 

III. PHYSICAL PROCESSES 

We focus on relativistic coUisionless shocks. Magnetic 
energy dissipation, for example via magnetic reconnec- 
tion, can also accelerate the particles and then give rise 
to prompt m, m, m, \^ and afterglow emission |47l. lisj . 
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FIG. 3: The prompt keV— GeV 7— ray emission (i.e., phase-1) and the long lasting high energy afterglow emission (i.e., phase-2) 
of GRB 940217 (from [s^). The source was earth-occulted for ~ 3700 s and the spectrum of all photons is inserted. 



However in this case we expect that the inverse Comp- 
ton emission will be weak because of the strong magnetic 
field (see eg. ip^ below) and it would not give rise to a 
significant high energy emission. 

In all scenarios that are considered in this review, the 
emitting material has been accelerated and it moves rela- 
tivistically relative to central engine with a bulk Lorentz 
factor r. Consequently there are two inertial frames: the 
rest frame of the emitting region (the comoving frame) , 
and the rest frame of the central engine. The observer's 
rest frame is the same as the latter aside from a cosmo- 
logical redshift factor (see also [13]). We denote param- 
eters measured in the comoving frame with the super- 
script ""' . The physical quantities (e.g. scale length and 
time) as viewed in the two inertial frames are related via 
the Lorentz transformations. Along the emitting region's 
moving direction in the central engine's frame, the length 
scales A' and A, and the time intervals elapsed for the 



same pair of events dt' and dt are related by 

A' ^ TA, dt' ^ di/V. 

t is the time measured in the central engine frame. This 
time is not that useful as we will be generally interested 
in the observer's time, that is the arrival time of photons 
as seen by an observer at infinity. For matter moving 
directly towards the observer the arrival time of a photon 
emitted at R relative to a photon emitted at i? = is: 



^obs — 



dr 



(l-/3r) 
c 



R 
2f2^ 



(1) 



where (ir = (1 - l/r^)"^/^. A general form (l-/3rCOs6') 
should be used in the integral above if the motion is at 
an angle 9 towards the observer. Since all observations 
are done with iobs we will use hereafter this observer time 
and drop the subscript "obs" sot = tohs- A cosmological 
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FIG. 4: The keV-GeV 7— ray light curve (left panel) and spectrum (right panel) of GRB 941017 in five time intervals from —18 
s to 211 s (from ;23]). 



facto (1+z) should be added to the above relation, as 
appropriate. Note that we have (for /3r = const.) 



i/2r^ = t'/r = t. 



(2) 



7ud(l + cr), fWT'A — \ — 2^ 
^ (7 + X + V 1 + 14x + X^) 

[1 ] 



(4) 



#-4> (5) 

A. Relativistic shocks 



Shocks involve sharp jumps in the physical conditions. 
Conservation of mass, energy, and momentum determine 
the Hugoniot shock jump conditions across the relativis- 
tic shocks. For the perpendicular shocks, if the upstream 
matter is cold and magnetized, the jump conditions read 



7 + X+\/l + 14x + X^ 



K 1 + X + + i4x + x^ 



(3) 



where n'^^, e'^^ and B'^ ^ are the number density, the 
energy density and the magnetic field measured in the 
local rest framed of the upstream (region u) and down- 
stream (region d), 7ud is the Lorentz factor of region u 

,2 

relative to region d, x= , < /c < 1 is the parameter 
describing the magnetic energy dissipation at the shock 
front (fc = 1 for the ideal MHD limit) and a is the ra- 
tio of the magnetic energy density to the particle energy 
density measured in the rest frame of the region u. 

For the un-magnetized upstream (i.e., cr = 0), we have 
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4 « 47.d, 
n, 



7ud 



1. 



(6) 



The fraction of thermal energy and the magnetic field be- 
hind the shock are needed to calculate the radiation spec- 
trum. These parameters are determined by the micro- 
scopic physical processes and are hard to estimate from 
first principles [8|. A phenonienological approach is to 
define two dimensionless parameters, Eb and Eg, that in- 
corporate our ignorance and uncertainties [5lj. So the 
energy of the shock given to the electrons and electrons 
are 



and 



respectively. 



U'^ = 4ec7ud(7ud - l)n„TOpC^ 
C/g = 4eB7ud(7ud - l)n„TOpC^ 



B. Particle acceleration 

Particle acceleration can occur in GRB blast waves 
through a First-order Fermi mechanism involving inter- 
nal or external shocks, and through second-order Fermi 
acceleration involving gyroresonant scattering of parti- 
cles by magnetic turbulence in the magnetic field of the 
blast wave. In the First-order Fermi acceleration, the 
particles are accelerated when they repeatedly cross a 
shock. Magnetic-field irregularities keep scattering the 
particles back so that they keep crossing the same shock. 
Using the relativistic shock jump conditions and kine- 
matic considerations one can find that the energy gain 
in the first shock crossing is of the order 52|. How- 
ever, subsequent shock crossings are not as efficient and 
the energy gain is just of order 2 ^5^. Repeated cycles 
of this type (in each of which the particles gain a fac- 
tor of ~ 2 in energy) lead to a power-law spectrum with 
p ~ 2.3 for 7ud S> 1 [11] ■ The Second-order Fermi accel- 
eration can work efficiently if magnetic turbulence in the 
magnetic field of the blast wave has been well developed 
[H^. Again, the accelerated spectrum could be a power- 
law with an index p ~ 2.2 [56]. In both scenarios, the 
maximal energy of the protons, limited by the strength 
of the shocks and the synchrotron cooling, might be as 
high as - 1020 eV [nisi [53. 

As mentioned above, in the shock front, Ec frac- 
tion of shock energy has been given to the fresh elec- 
trons swept by the blast wave. If the fresh elec- 
trons have a power-law energy distribution dn' /dj'^ oc 
{Ye ^ 1)^^ 7m < 7e < Ym: ^^^^ the shock jump 
conditions that J^^'' {dn' / d'y'^)dY — n'^ and 
l)mcC^{dn' /dYe)dYe = U'^, we have 



l-n 



3(7ud - 1) 



P- 



2 m-n 



1 me 



1, 



(7) 



where the maximal Lorentz factor is limited by the syn- 
chrotron losses and is given by [sll 

I'm « i^Y" « 4 X lO'B'-"\ (8) 

where e and are the charge and the Thompson cross 
section of electrons, respectively. The magnetic field 
strength 

B' = ^JsttUq = y^327r£B7ud(7ud - Ijn'^nipC^ 

- 0-04(^)'^'<'^'7ud/3ud Gauss, (9) 



where /3ud = ^1 — 1 /7^ ^ . So in the internal shocks or 
the early forward shock, B' is larger or even significantly 
larger than 1 Gauss as n'^ ^ 1 in internal shocks and 
7ud ~ r ^ 1 in the early forward shock. 



C. Radiation processes 

1. Synchrotron radiation 

The typical energy of synchrotron photons as well as 
the synchrotron cooling time depends on the Lorentz fac- 
tor of the relativistic electron under consideration and on 
the strength of the magnetic field. If the emitting region 
moves with a Lorentz factor F, the photons are blue- 
shifted. The typical photon energy in the observer frame 
is given by [5i] 



eB' 



^syn 



l + z 2t: (l + z)mcC 



(10) 



where 1?= [r(l - cos6')]-i. 

The power emitted, in the comoving frame, by a single 
electron due to synchrotron radiation is 59| 



-Psyn = ^CrTcC/B(7f - 1). 



(11) 



The synchrotron cooling time of an electron with a 
Lorentz factor 7^ is t'^ ~ "f'e'^^cC^ I Psyn- If is smaller 
than the dynamical time t'^ k, R/Tc, the electron cools 
rapidly (i.e., fast cooling), where R is the radius of the 
shock front to the central engine. We thus can define the 
cooling Lorentz factor of the shocked electrons (7e, which 
satisfies t'c — t'd) as 



7c 



"oIrb^ 



If 7e > j'm ' 11^*5 cooling of most electrons are not very 
fast, we call this case as slow cooling. The synchrotron 
radiation of electrons in turn modify their isotropic- 
equivalent energy distribution . The continuity equa- 
tion of electrons reads 



dNy^ 

dt' 



97^ ^ dV ■ 



(12) 
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where # « -^7^ and Q c j^' for < < 
7'm- 

dN , 

For a stationary distribution, -gp^ = 0. Wc thus 
have g^(^7^^) = Q- In both the fast and slow cool- 
ing cases, the electrons having a 7^ > 7c get cooled 
rapidly, for which the distribution should be iV^' cx 

{dYJdtr'lQdY « 7r^'^''- For < ie < ic, 
the cooling is unimportant so the distribution is still 
oc Q oc 7e~^. For 7c < 7e < 7m; we have (5 = and 
then A'-y/ oc 7e~^- Overall, we have the following distri- 
butions: 



7^ ^^^^\ for 7^ > max{7^, 7,'„}, 
/-p w ^/ \ ^/ \ ^/ (13) 
2 



iP\ for7^>7^>7;„. 



7e 



for 7m > 7e > 7c- 



The synchrotron radiation spectrum can be easily es- 
timated. The spectrum of one electron moving in a mag- 
netic field B' can be approximated by 



F'ix) « 2.149 



/3. 



(14) 



where x = v' jv'^^^ and i/' = (1 -f zyv^^v (y is the ob- 
server's frequency). We see that F'{x) peaks at a; = 1/3. 
If the synchrotron self- absorption is unimportant, for 
V < min{i/c, I'm}, the emission is the sum of the contribu- 
tions of the tails of all the electrons' emissions F^, oc v^^^, 
where i^c = i^synijc) and i^rn = i'syn(7m)- higher en- 
ergy range, using Fi.di' (x Nji_Ply^dy^ and i/ oc 7^^, we 
have Fiy oc for z^c < J' < ^^m, F^y oc for 

I'm < V < i^c and -Fjy oc v^"^!"^ for maxjz^m, ^^c} < i'- 
In summary, the synchrotron radiation spectrum can be 
approximated as (see also [60j ) 

for vm> V > maxjz^c, ^m}, 

j,-(p-l)/2^ for Vc>V> I'm, 
I'-i/^, for i/m > I' > fc, 
i^i/^, for min{zyc, i^m} > v. 

(15) 

The maximal specific fiux is estimated as j60l [6l| 

(1 + z)e37Vc,totrB' 



1^,3711— max 



where iVg^tot is the total number of electrons and is 
the luminosity distance of the emitting source. 

The maximal synchrotron frequency can be estimated 
by comparing the synchrotron cooling time with the ac- 
celeration time [SSl 



7Mr 



1 + z 2t: {I + z)mcC 



1 + z 



16(H-z)7r2e2 

In the prompt emission phase, F ^ 100, we have hvyi ^ 
a few GeV. In the early afterglow phase, /ii/m ~ 10 - 
100 MeV. If more energetic photons have been observed 
other mechanism(s) should be present. 



Protons can also produce synchrotron emission but this 
is, of course, much weaker as for 7p = 7^, the synchrotron 
radiation power of a proton is weaken by a factor of 
{iric/m-pY of an electron. We assume an initial power-law 
distribution dn' /d'^'^ oc (7p — 1)"'', for protons acceler- 
ated in the blast wave, the minimum Lorentz factor, the 
cooling Lorentz factor and the Maximum Lorentz factor 



/ m.p 



(1 



eB)(7ud-l)(p-2)/(p-l)-Hl, 



67rFmpC 



,RB 



/2 



rrir. 



7m, p 



(mp/me)7M, 



respectively, where the Thompson cross section of pro- 
tons is u^^-p = {mc/m-pYa^. For 7^ p > 7m p, we take 
7c p = 7J4 p. The synchrotron radiation frequency of the 
protons is 



^^syn,p 

l + z 



eB' 



27r(l -|- z)m-pC 



7p r. 



With Vr, 



^m,p — Z'syn,p(7m,p), i^c,p = '^syn.p (7c,p ) and I^M^p = 

i'sj/n,p(7M p)' straightforward to obtain the spectra, 
which take the form of cq. fTS]) . We note that 

VM,^, ^ {mp/mc)vM 

could be in TeV [g^]. But for reasonable shock param- 
eters that fit the afterglow data, the synchrotron TeV 
emission of protons is not as important as that of the 
synchrotron-self Compton of electrons [6^. In this re- 
view, we do not discuss it any more. 



2. Inverse Compton scattering 

An electron moving relative to a dense soft photon 
background will lose some of its energy via inverse Comp- 
ton scattering 59] and produce an inverse Compton com- 
ponent at higher energies 



2r j'^i^L 



1 + z 1 + z 1 



(17) 



where v!,^ is the frequency of the seed photon and g 
In the Thompson regime, g ^ 1, so 



i^ic ~ 7^ i^se- (18) 

In the Klein- Nishina regime, 5 > 1, we have 

« 7>ecV/i- (19) 

In this case, apart from the reduction in energy boost, 
the cross-section for scattering is also reduced to [s^ 
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3 .(1 + 5) 



2g(i + g) 

(1 + 25) 



ln(l + 2.g) 



(20) 



For convince, we define A{g) = a{v'^^,^'^/aT. 

The effect of inverse Compton scattering depends on 
the parameter 



y. — p! I pi 



syn 7 



(21) 



where is the power of the inverse Compton radiation, 
which can be estimated as « A{g)(TTc{Y/-l)Ui^/{l + 
g), where U!y is the energy density of the seed photons. 
We then have [6J| 



important. Note that second order IC will be even more 
important (see ea. (j36p below) and so will even higher 
orders. This divergence will be stopped by the Klien- 
Nishina cutoff. 

In general 



r°° HN' 



(23) 



Ajg) U', 
l + gUl,~ U 



1, 



for .9 < 1, 
for g > 1. 



(22) 



If yjc < 1, the inverse Compton effect is unimportant 
and can be ignored. On the other hand if Yic > 1 IC is 



The quantity dN'^/dt'du'.^ is the scattered photon spec- 
trum per electron [65j . Supposing the seed photons are 
isotropic in the rest frame of the IC scattering region, we 
can express dN'^/dt'dv[^ as (see ea. ([57|) for anisotropic 
photons): 



J 



dN'^ 
dt'dv' 



47f ^so 



[2qlng+(l + 2g)(l 



' 2 1+4^^ 



-?)], 



(24) 



where / = hv',^/ {-^'^meC^) satisfying hv'^^/ {"f'^rrieC^) < f < 
45/(1 + 45), q = //[45(1 - /)], and n'^,,^ is the frequency 
distribution of the seed photons in unit volume [H, [HHl . 

The cooling of electrons caused by synchrotron and IC 
radiation and adiabatic cooling is described by 

§=-:^-^[H-llchf-|, (25) 
dn 3mcC^ prl n 

where dR « VPrcdt'. Correspondingly, the cooling 
Lorentz factor takes a new form 



7c 



GTrFrripC^ 



2 ' 



(26) 



which is used throughout the rest of the review. 

Considering the spherical curvature of the emitting re- 
gion, the observed IC emission flux is 



(1 + ^) 

(i + ^)r 

4nDl 



dt'dv' ^' 



, dN' 
dtdv' 



d-f'jn 
yM^ (27) 



where is the solid angle. 

For electrons having a power-law energy distribution 
iV^' (X 7g ^, the IC spectrum is only weakly dependent 



on n' 



and can be approximated as (see eq.(2.76) and 
eq.(2.88) of [65^ for details) 



'(p-l)/2 



for g < 1, 
for g > 1. 



(28) 



Following [2J|, we present two different approaches to 
calculate the IC scattering with Klein-Nishina suppres- 
sion self-consistently. 

Instantaneous approximation. In this approach we as- 
sume a functional form for the energy distribution of the 
electrons acceleration in the shock front, n'(7g), and con- 
sider its instantaneous modification due to cooling. An 
electron of Lorentz factor 7^ has a cooling time given by 



(29) 



If ic(7e) ^ ^d; then the electron emits both synchrotron 
and IC radiation for the entire time t'^. However, when 
*c(7e) < the electron radiates only for a time tc(7e)- 
Thus, the total spectral radiation density, including that 
of the seed photons and that produced by all the electrons 
in the fluid, is given by: 



+ P4(^',7e)] xMin[i:„i^(7^)]n'(7^)d7^, (30) 
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where Ply^{iy',ji) « {aTm,c-'B')F'{i,'/iy' )/3e and 
^;(^',7D«(l + 5Wa(^',7D- 



Equation (|30p is an integral equation, since the func- 
tion Pl^{v' ,"f'^) inside the integral itself depends on [/^,. 
The quantity 7,'„i„ is the smallest 7^ down to which elec- 
trons are present. In dealing with equation pop we need 
to consider two cases. One is the slow cooling^ in which we 
may use equation ([50)1 directly with 7^;^^ = 7^ and f^'(7e) 
given by the original energy distribution produced in the 
shock. The other is the fast cooling, in which electrons 
will continue to cool below 7,„ to a minimum 7|^ijj such 
that t^(7min) = irf- Now, for the range 7^„ < 7e < 7™, 
all the electrons are available for radiating. Initially, most 
of the electrons are at 7^, and as these electrons cool each 
electron will pass every 7^ between 7^ and ^'^^^^ (where 
all these electrons accumulate). Hence we have 



ri'(7e)-«'(7;j, 



Till in 



< 7e < 7™ 



(31) 



As usual, we assume a power-law distribution for the 
electron Lorentz factor n'{Y^)djl cx j'e~^dj'^ (7^ > 7^), 
for which 7^ is given by eq.([7]). Equation pop may be 
solved numerically via an iterative method. The algo- 
rithm proceeds as follows. We begin with some reason- 
able initial approximation for Ul,. Using this, we com- 
pute i^'j,(7e), tc(7e) ^iid 7^in. Then, we compute the 
spectral distributions Pgy„(z^',7g) and P(^{i^' , ■y'^) for all 
7e ^ 7min obtain via equation ([50)1 a new approxi- 
mation for Ul,. We take this U^,, or (for smoother con- 
vergence) a suitable linear combination of the new and 
old U[,,, as the current approximation for Ul, and repeat 
the steps. The iteration usually converges fairly quickly. 

Dynamical approach. In the limit of a single IC 
scattering, with given n'^, , T, B' and Q, eg. p^ . 
ea. ipT]) and eqs. (|23ll2"5)) form a complete set of formulae 
which can be solved for N^i numerically. Then the 
IC flux can be obtained using eg. ([27)1 . The treatment 
is a bit more complicated if the synchrotron radiation 

component n^, « -^^j^^F' {v'J v'^y^We (see 

eq. (fT4ll for i^'(i/'/i/gyjj)) is comparable to the seed photon 
background. In this case we need to replace the term 
n'^, in eg. ((24)) by rt^, +n'^, \u'^=v^^ and then solve for N.yi^ 



self-consistently. The resulting n'^, 



J dt'du' AttR^c 'e 



dt' du' 4-71 

in turn plays an role in cooling the electrons. Combining 
this IC component with n[^, and the synchrotron 
radiation component of electrons, we can solve for N-^i 
self-consistently and then calculate the second order IC 
spectrum and so on. 



parameter Y^sc can be estimated as 



Y 



A{g) Ujy, _ 

1+g U'g l + .g''£B l + i"s 



A{g) e. 



(32) 



where 77 = min{l, (7^/7^)^^^} and the relation Ui.y^ — 
riU^/{l + Yssc) has been taken into account. In the 
Thompson regime, we have [67j 



Y 



-1 + ^1 + 4?/£c/£B 

2 

■qSc/eB, for r]ec/eB < 1, ,„„^ 
(?7ee/eB)l/^ for ve,/eB » 1. ^ ' 



The SSC spectrum can be approximated as follows: 



for I^M.ssc > 1/ > max{z^c,ssc, ^'m,ssc}, 



. V for I/^ssc > V> Vm,s 



> V > Vr. 



,1/3 



for v„ 

for V < min{i^, 



c.ssci ^m.ssc 



where Vm.ssc ~ '^I'mV^, Vc^ssc ~ '^I'cVc [68|, and the SSC 



,,'2, 



(34) 



cut-off frequency vm,s 



above which the 



IC is in Klein-Nishina regime and is very weak. With 
a given synchrotron spectrum -Fjy.syn taking the form of 
eg. ifTS)) . the SSC spectrum is fixed by 



Y 



(35) 

If max{7^ Vm, 7^ < rn^.f? (2, the second order 

IC scattering is still in Thompson regime. The luminosity 
ratio is given by 



Y — 

^ SSC — 



1 



U'b 



l + ^ss 



2ndIC Ug 



where F^ndic = U'^JU'b = {UUUly^){UlyjU{,) = Y^^, 
which dominates Fic if Ksc > 1- For ijU^/Uq = tjEc/s-b ^ 
1, we have (see also [sl, [t^I) 



Y 

^ ss 



(?7eo/eB)^^^ >2ndic ~ [riec/eBf^^- 



(36) 



Most of the energy of electrons are lost in the second scat- 
tering. Now t'm,2ndIC ~ 47^^f„ and I^c,2ndIC ~ 'i'lc'^'^c- 



The spectrum takes the same form of that of eq. ([34|) . 
Similarly if the third scattering is still in Thompson 
regime and rjEc/eB ^ 1, we have Fssc ~ (?7£c/£b)^^^ and 

>^3rdIC ~ {neJSBf'^. 



Synchrotron- self Compton (SSC) takes place if the 
seed photons involved in the inverse Compton scattering 
are the synchrotron radiation of the shock-accelerated 
electrons (i.e., U'^ — U!.y^J. In the comoving frame of 
the shocked material, both the electrons and the seed 
photons are isotropic. The Synchrotron-self Compton 



External inverse Compton (EIC) takes place if the 
seed photons are from a region well separated from the 
scattering electrons and along the direction in which the 
ejecta moves. In this case, the electrons are isotropic 
but the seed photons, of course, are highly beamed. The 
spectrum of radiation scattered at an angle ^sc relative 
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FIG. 5: The normalized emitted power as a function of the 
scattering angle for different power-law energy distribution 
electrons (from [t^])- From bottom to top, the power-law 
indexes are 4, 2.5, 2.0 respectively. 



to the direction of the photon beam penetrating through 
this region is 7lj : 



dt'dv'. dVL' 



3aTC '^'u'^^dvi.^ 
2^ . 



e 



2(1-0 



2e 



he{l-i) hl{l-if 



(37) 



where dO,' = 27rsin0scrf^sc, £, = hv' . J {'y'^mf.c'^) , bg = 
2(1 - cos9scWehvU{meC^), cos6lsc ="(cos6' - /3r)/(l - 
/3rcos0), and hv!,^ <^ hv'^.^ < "f'^rrieC^he / {1 + bg). In the 
case of isotropically distributed photons, the averaging of 
eg. (157)1 over the angle ^sc reduces to ea. ip^ . as expected. 

In the rest frame of the emitting region, the EIC 
emission has a maximum at ^sc = tt and it vanishes for 
small scattering angles, as shown in Fig [5] This effect 
lowers the EIC flux in two ways 2j|. First, a fraction 
of the total energy is emitted out of our line of sight 
and thus the received power is depressed (relative to the 
isotropic seed photon case). Second and more important 
for GRB study, the strongest emission is from 9 ^ l/T 
(see eg. ([55]) below). As a result, the EIC emission 
duration is extended and the flux is lowered. This will 
have important implications on EIC from X-ray flares as 
discussed in section HV A 31 below. 

Bulk Compton Scattering. Inverse Compton may 
arise even if the electrons in the ejecta are cold but are 
moving relativistically into a soft photon background. In 
this case, the IC spectrum depends on both the initial 



spectrum of seed photons and the deceleration of the 
electrons. Roughly speaking, the upscattered photons 
have a much larger frequency v^ic ^ ^^v^c and the flux 
is Fy^.^ ~ Tbic-f^i/.e , where Tbic ~ crTA^e,tot/(47ri?2), v^c is 
the frequency of the background photons (measured by 
the observer) and F^^^ is the spectrum. 



3. J— rays from pion production 

High energy photons can also be produced via 7r° decay 

TT*^ ^ 7 + 7. (38) 

In GRB internal shocks and blast waves, 7r° can be pro- 
duced via the following processes: 

p+p — > p+p+tt'^ , p+n — !■ p+n+n'^ , P+7 ^ ir^+p. 

(39) 

Please note that in this sub-subsection, p and n represent 
proton and neutron respectively and 7 represents photon. 
In the rest frame of the outflow, the inelastic pion pro- 
duction threshold is e 



pion 



140 MeV. For the observer, 
the resulting 7— rays thus have an energy > TOT MeV. 
Actually, the 7— rays resulting in p+j A+ — > ir'^+p — > 
7+7+p are much more energetic because the cross section 
of p -I- 7 ^ A+ peaks when x eproton ~ (0.3 GeV)^r^, 
where is the typical energy of prompt 7— rays and 
Eproton IS the energy of VHE protons. So wc have 

eproton ~ 0.1 rV[e7/(l McV)] TeV. (40) 

For the observer, the photons from the 7r°— decay is as 
energetic as ~ 5 r^/[e-,/(l MeV)]GeV if the energy of 
tt'^ is ^ 10% of the incident protons. 

Alternatively, high energy photons are also arise from 
in the synchrotron radiation of the pairs that result in 
the , TT^ decay 

r + 



IT 

7r~ 



9^ 



(41) 



Possible processes involved in GRBs to produce 



TT"*", TT are following 



P 



p + n 

p-f 7 -H 



p + p^p + n + TT'^ 
\- p + . 
TT ' + n. 



A+^7r++n. (42) 

The last process of the above four is particularly inter- 
esting because the resulting e"*" has a very high ran- 
dom Lorentz factor ~ 0.05 eproton/(rr7ieC^) ~ 
10^r/[e-y/(l MeV)]. It's synchrotron radiation may peak 
at TeV energies. Though interesting, the energy of this 
TeV emission component is not expected to be more 
than that of the synchrotron soft 7— ray emission un- 
less Ec < 0.01 (e.g., [HI). This is because the fraction 
of total shock energy given to protons above eproton ~ 
10^ rl 5/[e.^/(l MeV)] GeV is only < 1/3. Furthermore, 
only a small fraction (< 0.1) of the energy of these ultra- 
relativistic protons is lost in producing e+. 
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4- Electromagnetic cascade of TeV 'y—rays 

As the 7— rays with an energy ~ ITeV travel toward 
the observer, a significant fraction of them will be ab- 
sorbed due to the interactions with the diffuse infrared 
background, yielding pairs [H, H^. If a primary pho- 
ton with energy has been absorbed, the resulting 
pairs have Lorentz factors 



7e 



y/(2meC^) 



10^ e. 



/(ITeV). 



(43) 



Such ultra-relativistic pairs will subsequently (bulk) 
Compton scatter on the ambient cosmic microwave back- 
ground (CMB) photons, and boost them to an average 
value 



0.63(1 + z)(e^/l TeV)2 GeV, (44) 



where e = 2.7fcTcinb is the mean energy of the CMB pho- 
tons with Tcmb ^ 2.73(1 + z}K. 

As shown in [tI, [tI, [t^, 0] , there are four timescales 
involved in the emission process (all are measured by the 
observer): The first is tact, the observed activity time 
of the source emission. For GRBs, it is unlikely to be 
longer than 10'^ sec. The second is the well-known angu- 
lar spreading time 



At A 



(1 



= 960(1 + z)(-^)-^(—:i^^ 
27^c 10° 0.1cm 



where Rr. 



(0.26crTnj, 



5.8xl025(_^ 



is the typical pair-production distance, and n,-^-^ ~ 
O.lcm"^ is the intergalactic infrared photon number den- 
sity. The third is the inverse Compton cooling timescale 



At 



ic 



(1 + z)tic,loc/(27e) = 38(1 + z)-3(^^/106)-3 g. 



where the IC cooling time scale measured in the source 
frame iic.ioc = 3mec/(47e(TTMcmb) = 7.7 x 10^^(1 -I- 
z)-4(7,/106)-i s, Ucmb = a7;U is the CMB energy den- 
sity, and a is the radiation constant. The fourth is the 
magnetic deflection time that arises due to the deflection 
of the pairs by the intergalactic magnetic fields (IGMF): 



6.1X 10"(^^-^ 



10^' 4o-i6G 



where Biq is the strength of the field. Clearly this delay 
is much too long unless Biq is extremely small. 

In the most optimistic case, the prompt or very early 
GRB afterglow flucncc in TeV energies would be 5tcV ^ 
10~^ erg cm^'^, thus the expected GeV flux is (assuming 
AtB dominates Atdeiay) 



At 



TcV 



delay 



3 X 10-13 (^)5( 

Ho^' 4o 



B 



IG 



^16 



G 



TcV 



10 



■^erg cm ^ 



) erg cm ^ s 



(45) 



-SiG < 10^1® G is required to give rise to detectable sig- 
nals (see also [t^, [zll). However, a recent estimate of 
Big gives Big 10"^^ G [7§| and values as large as 
0.1/iG have been suggested. For such an IGMF, the cas- 
cade radiation from the pairs will appear as a halo 
around the TeV 7— ray source because the magnetic fleld 
is strong enough to make the distribution of these pairs 
isotropic [tI]. As a result, the flux will be much too low 
to be detectable. 



IV. HIGH ENERGY EMISSION PROCESSES IN 
GRBS AND AFTERGLOWS 

We begin with a brief discussion of the standard the- 
oretical model of GRBs and afterglows (see FiglH] for a 
schematic plot) and refer the readers to Q for a detailed 
review. 

• The bursts are produced by compact sources 
that generate relativistic outflows — relativis- 
tic jets with initial bulk Lorentz factors F 

O 



tens - hundreds [8l|, US, HI . 

• The emission is produced by shocks that accelerate 
particles and generate magnetic fields. An alterna- 
tive interpretation is that magnetic dissipation via 
plasma instabilities and reconnection tak es p lace 
in a Foynting-flux dominated outflow [H, [ij . 
The shocks (or the magnetic dissipation process) 
can be internal (if they are within the outflow) 
[H, m, m, 113] or external (if they are due to in- 
teraction with the surrounding matter). Note that 
internal shocks must take place at a radius smaller 
than the deceleration radius (^ 10^^ — lO-'^^cm), in 
which the outflow is slowed down by its interaction 
with the surrounding matter. Note that also in the 
case of magnetic dissipation case timing arguments 
suggest that the dissipation is due to internal pro- 
cesses. 

• When the outflow is decelerated by the surround- 
ing medium two distinct shock waves form [88| . 
A forward shock that expands into the medium. 
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FIG. 6: A schematic plot of internal-external shock model for 
Gamma-ray Bursts and their afterglows (from [80l|'). 



prompt 




E 

VII'. I 



XRT lightcurve 



log (t) 



and a reverse shock that penetrates into the out- 
flow. The typical Lorentz factors of the (forward, 
reverse) shock electrons are ^ (10'*, 10^), respec- 
tively. We denote the phase during which both for- 
ward and reverse shocks exist as the very early af- 
terglow. In this phase, the forward shock emission 
peaks in X-ray band while the reverse shock emis- 
sion peaks in the infrared or optical bands [89] . The 
reverse shock may lead to an optical flash that be- 
gins shortly after the onset of the prompt 7— ray 
emission. Such a flash has been detected in several 
powerful bursts l91|]. After the reverse shock 
crosses the outflow, it dies out and only the forward 
shock remains. A self similar blast wave forms that 
propagates into the surrounding material. For a ho- 
mogenous circumburst medium, the bulk Lorentz 
factor of the blast wave is ISOl 



6El%n-'/\t/l day)-3/«(l + z)^/«. 



(46) 



where -Ek is the isotropic-equivalent kinetic energy 
of the ejecta, and n is the number density of the 
medium. Throughout this review, the convention 
Gx = in cgs units has been adopted. The 

synchrotron radiation of the shock-accelerated elec- 
trons [51I, [9^. .931 fits the late-time GRB afterglow 
pretty welllSflH ^ , though at times jets , 
a wind medium 98, lOCrl, and energy injection 
fiol I1O2, 103, 104, 105] had to be invoked to ac- 
count for the observations. 

• Prolonged activity of the GRB central engine 
also plays an important role in producing after- 
glow emission (i.e., the central engine afterglow) 
llM IToi either via late internal shocks ^M, flTol 
lllll IllsT 1 1 3[ or via late magnetic energy dissipa- 



tion [^, [4^, III4 1 . This activity has been suggested 
to interpret the peculiar features that emerged in 
the Swift's observations of the early {t < 10^ sec) 
X-ray afterglo w (see Figl?]): energetic X-ray flares 
[lllllll5llll6ll. X-ray plateaus that are followed by 
sharp drops [83, Ill7l | and X-ray flattenings asso- 



FIG. 7: A schema tic p lot of the X-ray afterglow of GRBs (see 
also Zhang et al. [l06l | and Nousek et al. flO^ ). 



dated with ch romatic breaks in the optical band 
[TT8l . [lT9lfl20l . 



A. Inverse Compton 

Inverse Compton i n the for m of either SSC [gl, Il2ll 
m, [123 or EIC m [ill [m is Hkely the most impor- 
tant source in producing high energy 7— ray emission. 



1. Internal shocks SSC 

For internal shocks taking place at a typical radius R, 
the magnetic field can be estimated as 



B' ^ [2{l + Y,,,)eL,yJ{rlR^c)] 



]l/2 

-1/2 



~ 10^ Gauss (1 + Y,,,)-'^'Ll;^l,,T-'R^,\ (47) 

where e = ee/ee ~ (1 + i^sc)~^, and Lgyn is the syn- 
chrotron radiation luminosity of the internal shock emis- 
sion. The corresponding typical electron Lorentz factor 
is 



7e 



UOO{l+Y^.,y^'L:'J'Rli'a+^)'^'i 



• 100 keV 



1/2 



(48) 

where Cp = hv^^ is the observed peak energy of the syn- 
chrotron emission. The energy of a typical inverse Comp- 
ton photon is 



hv„ 



^ 240GcV (1 + Yssc) 



1/2 



100 keV 



(49) 



We summarize the typical values of the parameters in- 
volved in ea. (|49)) and the expected peak energy of the 
SSC emission in Table Hill The SSC spectrum of internal 
shocks of brigh t GRBs can show a significant GeV— TeV 
signal psi . Il26j . as shown in Fig [5] 
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TABLE III: SSC emission of internal shocks in the prompt emission phase a nd in centr al engine afterglows. Lssc is the SSC 
radiation luminosity of internal shocks. Note that Lssc/Lsyn ~ ^sc ^ ( — 1 + a/1 + 4ec/eB)/2 ^ 1 for Ec/eb ~ a few. In the last 
column Y, P, N represent Yes, Possible and No respectively. 



phases 


-^syn.SO 


Rvs 


Ep/keV 




-^ssc /-^syn 


Detectability 














LAT (MAGIC) 


GRB 


1 - 100 


0.1 - 10 


10^ - 10^ 


GeV - TeV 


~ 1 


Y (P for z <1) 


X-ray Flash (XRF) 


10"^ - 1 


0.1 - 10 


~10 


GeV 


~ 1 


Y (N) 


Central engine afterglow 


10-5 - 0.1 


10 - 10^ 


~0.2 


sub - GeV 


~ 1 


P (N) 




Obs. energy [eV] 



FIG. 8: The GeV excess due to the SSC emission in internal 
shocks: the synchrotron + SSC spectra of the GRB internal 
shocks (from [25l]). The synchrotron + SSC luminosity L — 
10^2 erg/s, Ee = 1/3, r = 600 and SU = 10"^ s for aU lines. 
From top to bottom (i.e., the solid line, the dashed line and 
the dotted line), eb ~ (1/3, 10"^, 10"^), respectively. 

Prompt VHE photons above the cut-ofF frequency Vcut 
will produce pairs by interacting with softer photon s and 
will not escape from the fireball. Following [tIi Il27j | , 
we have 

hi^.ut - 2 GeV (ep/100 keV)(2-rt/Pi;^2/P^ 

stl%ri'',l''^r (50) 

where is the observed variability timescale of the 
prompt soft 7— ray emission light curve. Consequently 
the prompt TeV emission can escape only if To > 10^. 
Additionally, as mentioned in section III Bl the Universe 
is optically thick for TeV photons. Hence TeV emission 
could be detected only for bright (high To) and nearby 
GRBs, which of course are quite rare. 

For X-ray flashes (XRF) and X-ray rich GRBs, the 
peak energy is so mewhat lower than in GRBs but still 
above GeV [l28j |. If the X-ray flare photons are due 
to the synchrotron radiation of late internal shocks, the 
SSC emission component will peak at GeV energies and 
it will be de t ecta ble by the upcoming GLAST satellite 
[H, [m HH mH . K x-ray flares are powered by the re- 
freshed shocks (produced by slowly moving matter that 



was ejected more or less simultaneously with the faster 
moving one during the onset of the prompt emission) 
at a radius ~ 10^^ cm, the SSC emiss ion m ay peak in 
GeV- TeV energy range, as shown in [2^.ll30l|. The X-ray 
plateaus that are followed by a sharp drop, such thos e de- 
tected in GRB 060607A [83] and GRB 070110 [uS. are 
likely to be the synchrotron radiation of shocks powered 
by the prolonged activity of the central engine. The SSC 
radiation components associated with this process ap- 
pears as GeV plateaus and we expect it to be detectable 
by LAT [133. 

2. Inverse Compton processes in the very early afterglow 

We denote by very early afterglow the short (~ 100 — 
1000 sec) phase during which both forward and reverse 
shocks exist. This gives a rich structure of possible inter- 
actions of photons from dil ferent regions with electrons 
from different regions [l23l . fisi [TsslflM [T37| . The sit- 
uation is more complicated if the prompt photons over- 
lap the reverse/forward shock fronts. Such an overlap 
is important if the reverse shock crossing radius satisfies 
-Rx ^ 5 X 10^^ cm Fq 2.5^0,11.5, where Aq is the width of 
the GRB ejecta [63|. The shock crossing radius i?x can 
be written in a general form 

i?x ~ max{i?.^, 2F2 Ao}, (51) 

where is the radius where the mass of the medium 
collected by the fireball is equal to I/Fq of the fireball 
mass and Fx is the Lore ntz factor of the shocked ejecta 
at the crossing time |l38{ . In eq. ((5T|) , the first term dom- 
inates if the reverse shock is weak (Fx ~ ro/2), while the 
second term dominates if the reverse shock is relativistic 
(Fx ^ To). Generally we expect that for a weak (rel- 
ativistic) reverse shock the overlap between the prompt 
emission and the external shock fronts is unimportant 
(very important). 

Inverse Compton emission can take place in the very 
early afterglow from either SSC or EIC. There are two 
cases depending on the strength of the reverse shock, (a) 
A weak reverse shock, in which the prompt soft 7— ray 
photons exceed the external shock fronts quickly. Its ef- 
fect on cooling the reverse/forward shock electrons is ig- 
norable. In this case, only a factor of ^ 0.1 of the total 
energy [133| | is given to the reverse shock and the rest 
is given to the forward shock. The seed photons for the 
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TABLE IV: Inverse Compton emission of very early afterglow 

Rev. shock electrons For. shock electrons Cases 

Rev. shock photons hum,ssc ~ 10 keV(r/300)'^ hum.,ic ~ 100 MeV(r/300)'' a, b 
For. shock photons /li^^^.ic - 100MeV(r/300)* /ii/m,ssc 1 TeV(r/300)® a, b 

GRB photons fei^m.cic ^ ^ G!eV(^^j^)^( -^p^^^y ) Klein-Nishina regime b 

XRF photons ^50MeV(4ig-)^(^) fci.„,,ic~5GeV(^)^(^) b 



reverse (forward) shock electrons are their synchrotron 
radiation and the forward (reverse) shock synchrotron 
radiation jl23l Il35j . With typical parameters the energy 
of the synchrotron photons and the electro ns' L orentz 
factor in the forward and reverse shocks are }l37| : 



J^syn.F ~ lOkeV (^/300)^ 

i^syn.R - 1 eV (r/300)2, 



7e,R 



104(r/300); 
100. 



As shown in Table IIVI there are four inverse Comp- 
ton processes that produce high energy emission. The 
most important high energy signature in this case is the 
GeV— TeV SSC emission of the forward shock. 

(b) A strong relativistic reverse shock, for which the 
overlapping between the prompt emission and the for- 
ward/reverse shock regions is tight. In this case, the 
reverse shock has an energy that is comparable to the 
forward shock an d the prompt emission overlaps the re- 
verse shock front [l38j |. The reverse shock electrons will 
be mainly cooled by the prompt ~ 100 keV 7— rays and 
give rise to EIC emission with a typical energy [l24| 



27' 



1 GeV(^)2( 



100' 'lOOkeV' 



In the forward shock region, 7^ pCp/F 3> nicC^ for 
Ep ^ 10^ — 10"^ keV and the EIC process is in the Klein- 
Nishina regime. So the forward shock electrons will lose 
their energy mainly via SSC unless the prompt X-ray 
component is so energetic that can cool the forward shock 
electrons effectively [64| . If the prompt emission is so soft 
that Cp ~ a few keV (i.e., XRF) and F ^ tens, the IC 
scattering of the forward shock electrons on the prompt 
emission is still in the Thompson regime and GeV EIC 
emis sion is expected. One good candidate is XRF 060218 
[l39l |. The results are summarized in Table HVl 

In the case of a dense wind medium, the reverse shock 
crosses the outflow at a radius Rx ~ 10"^^ cm <C 2FqAo. 
At such a small radius, because of the tight overlap of 
the prompt 7— rays and the forward shock, the optical 
depth for the GeV— TeV photons produced in the shocks 
may be very large (see eq.(12) of !64|). These high-energy 
photons interact with the prompt photons and generate 
relativistic pairs. These pairs re-scatter the soft X- 
ray photons from the forward shock as well as the prompt 



photons and power a detectable high-energy emission, a 
significant part of which is in the sub-GeV energy range. 
Consequently we will observe an energetic delayed (see 
eq.® below) sub-GeV flash [6^. 

However, bright optical flashes hav e be en detected 
only in quite a few bursts (e.g., [ool. [olI . [l40j ). for which 
the afterglow modelling suggests a wea kly m agnetized 
(cr < 0.1) reverse shock region ([Hi, [HO, [mI J 42, 143], 
see however [144 Il45l Il46| ). The non-detection of 
optical flashes in most GRB afterglows may imply a 
mildly or highly (cr > 0.1) magnetized outflow, in which 
the SSC emission is very weak. So the contribution of 
the reverse shock to the high energy emission may be 
unimportant in most cases (see however [TOjl. 



3. SSC and EIC in the afterglow 

During the late afterglow the main radiating region 
is the forward shock that moves into the surrounding 
matter. We consider both an ISM (k — 0) and a wind 
{k ~ 2) medium. 

SSC: standard afterglow model. The two charac- 
teristic frequencies governing the spectrum are 



10^1 Hz C*ei^ 



1/2 

l^B,-2 



6.2n-VX^53(l- 
1.4 A-i/i'i?k,53(l 



- z)t 



3 
-2 
3 ' 



for A: = 
for fc = 2 



(52) 



V.^ssc ~ 10^^ Hz (l+i;sc)"V-2 



4 n—E^'^^^{l + z) — t^'' , for fc = 0, 
1.5 A7:i£;k,53(l + zyHl for fc = 2, 



(53) 



respectively [631 [681, [1231 . Where Cp = 13(p - 2)/[3(p - 



1)1, A 



[M/lO-5M0yr-i]K/(lO8 



is the wind 



parameter, M is the mass l oss r ate of the progenitor, 

is the velocity of the wind |lOO( . 

The maximum SSC flux can be estimated as [6a Il47l 
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.n 11 2 1 . r w 1 [ 0.07 ^5/44/2 .El^^^tl^^m^f^DTl. 34, for fc = 0, 

I 1 /4^/ e ' f-^ ( 1±£ \2 far k - 9 

I ^ ^*,-l^B,-2''3 I 2 ^L,28.34' ™ K — Z. 
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FIG. 9: The multi-wavelength spectra of a GRB forward 
shock at 2 X (10^, lO", 10^) sec (from top to bottom). Thin 
and thick lines correspond to the pure synchrotron spectrum 
and SSC+synchrotron spectrum, respectively (from Be- 
cause of the Klein-Nishina correction, the VHE SSC spectrum 
is much softer than that of the synchrotron X-rays. 



With eas. ([34l) . we can estimate the SSC emission 
flux at a given frequency. The relation t] = 
min{l, (i^m/t^c)^^"^-'^^} is needed to estimate Fssc, where 



0.024(1 + z)Clel „i4 _2"^^k 53^3 \ for fc = 0, 
zYClel 



0.12(1 



''p = G,-l=-B,-2^*,-l''3 



for fc 2. 
(55) 

A numerical example of the high energy afterglow 
spectrum is shown in Fig [91 

SSC: modified afterglow model. 

The standard GRB fireball model fails to explain the 
shallow deca y phase of the Swi ft GRB X-ray afterglows. 
As shown in jnl [HI [Hi, [mI] , we need to consider vari- 
ous modifications: (1) The fireba ll underg oes a significant 
energy injection, i.e., Ek oc t^~'^ jl04ll05j |; (2) The s hock 
parameters are shock-strength dependent jllSl . Il48l |. i.e., 
£0 oc F"" and eb oc F"''. In both cases, an early flatten- 
ing is expected in the high energy afterglow light curve, 
as shown in [2^. [l30l [mtTITsI [Tsij . Without these mod- 
ifications, both Ek and £e take a constant value from 
the beginning. The early time X-ray emission would be 
stronger but drop with time faster. The weak X-ray sig- 
nal in the early afterglow phase implies a less optim istic 
detectio n pro spect of the GeV emission [13, [H, Il50l | (see 
however [l5ll |). 

Table |V] gives the spectral and temporal indexes 




strongest IC 
radiation is 
o from e ~1/r 



medium 



FIG. 10: A schematic plot of the external inverse Compton 
scattering: the case of the (late) prompt 7— rays. X-rays and 
ultraviolet photons from the central engine upscatter ed b y the 
shock-accelerated electrons in the blast wave (from [l52|| '). 



a and /3 of the afterglow emission. We define 
a = ao + ctE + + ay, where ao corresponds to 
the contribution of the standard emission, a e represents 
the contribution of the energy injection, stands for 
the contribution of evolving shock parameters, and ay 
comes from the evolution of Compton parameter l^sc- 
For example, if we only consider energy injection, then 
a = oq + OLE- If only evolving shock parameters is 
considered, a = ao + a^. If both effects are considered, 
a = ao + -|- a„. If Fssc ^ 1, the term ay should be 
included Il47l. 



EIC. It has been suggested that some GRB central en- 
gines produce late prompt emission in the form of X-ray 
flares or X-ray plateaus followed by sharp drops. These 
late prompt photons with an energy ex will catch up with 
the blast wave, cool the shock-accelerated electrons and 
give rise to EIC emission [125, ,152] , as shown in FigfTOl 
In this case, the peak energy of the EIC emission is 



hvn 



0.4 GeV 7t,3(ex/0.2keV), (56) 



where 



1.7 X 10^ Ec-iCp 



El[W%''\^ + ^)/2]^/«, for fc - 
z)/2]V4^ 



.1/4 4-l/4,-l/4r 



for fc = 2 



f57) 



A novel phenomena that appears in the EIC process is 
that the duration of the high energy emission is signifi- 
cantly longer than the duration of the seed photon pulse. 
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TAB LE V: The spectral and temporal index f3, a of afterglow emission in the case of ISM. Here (x v "is adopted (from 

[na). 
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FIG. 11: The EIC emission with (the dotted line) and without 
(the thin dashed line) anisotropic correction (from [IJ). The 
thick lines are the SSC emission of the forward shock in two 
different energy ranges. 



This is because the duration of the high energy emission 
is af fected by the spherical curvature of the blast wave 
[l24'| and is further extended by the high ly anisotropic 
radiation of the up-scattered photons |l52l | . One can see 
these two effects, particularly the second one, in FigfTT] 
Below, following f24!|, we give an analytical derivation 
of the peak time of the EIC emission caused by a seed 
photon pulse. For electrons with an energy distribution 
dn' /dy^ oc 7g if the IC scattering is in Thompson 
regime (i.e., ^ ^ 1), with eg. ([571) can show that in 
the local frame of the shocked medium, the resulting EIC 
emissivity is proportional to (1 — cos 6'sc)^^~''^''''^ [z3- The 



observed emission from an angle 9 is thus 

oc 0P+2(l+r202)-(7+p)/2, (58) 

The term T>^ is caused by Lorentz translation of the emit- 
ting power, sin0 is from the expression of the solid angle, 
and (1 — cos6'sc)'"'^^^''^^ is the anisotropic correction of the 
EIC. 

We define a 9c at which the emission peaks (i.e., 
^ \9=e,= 0) and have 

^c«(^)^/'r-i«i/r. (59) 

5 

The EIC emission thus peaks at a time 

rp«(l + z)i?(l-cos0J/c«^i±^, (60) 

where R is the radius of the shock front. 

For UV/X-ray flares having a duration 6t ^ 0.3if, the 
duration of the EIC emission is Tdc Tp ^ 3tf , where tf 
is the flare peak time. As a result, Toic ~ lOSt. 

The luminosity of the EIC emission can thus be esti- 
mated as 

cic -^p 

where L^^^ is the power given to the freshly shocked elec- 
trons in the blast wave and can be estimated as 

/ 7.5 X 10"* erg s-\ for fc 0, , . 
\ 5 X lO'^s erg s-\ for k = 2. ' 

We note that L^^^ depends only weakly on the density 
profile. 
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4- Bulk Compton 

Shemi [m and Shaviv & Dar [HI [HI proposed bulk 
Compton scattering as the source of the prompt 7— ray 
emission. In their scenario, the ultra-relativistic ejecta 
is moving into a dense soft photon background and the 
electrons in the ejecta Compton scatter on these photons 
and boost them to MeV energies (see also (l56j |V 

Bulk Compton scattering can also take place in inter- 
nal shocks, as suggested by Takagi & Kobayashi jl57l |. 
In this case, the seed photons are the synchrotron pho- 
tons emitted by one of the internal shocks. They are 
up-scattered by the electrons carried by faster shell(s) 
ejected at late times (before this shell collides so the elec- 
trons are cold). The bulk Compton scattering produces 
a lOOMeV-GeV emission component if the sub-outflows' 
Lorentz factor varies significantly, say, between 10 and 
10^. The efficiency of producing high energy emission 
in such a process, however, is low to ~ 0.001 — 0.01 for 
baryon-rich outflows. 

Recently, Panaitescu 'iSS] suggested that the X-ray 
flares, the X-ray plateau followed by a sharp drop, and 
the shallow decay X-ray phase that ends without an op- 
tical break are produced by bulk Compton scattering of 
the forward shock synchrotron photons by electrons car- 
ried in a new outflow launched by the central engine. In 
this model, the bulk Lorentz factor of the new outflow 
r„cw has to be much larger than the decelerating GRB 
outflow. In [Hi], Tnew is taken to be ~ 10'' - 10^ The 
new outflow has also to be pairs dominated otherwise 
the efficiency of energy conversion in the bulk Compton 
scattering process would be very low [l57j . It is not clear 
whether these two conditions are realistic or not. If cor- 
rect, a GeV emission component associat ed w ith the pe- 
culiar X-ray afterglow should be present |l59l |. Unfortu- 
nately, it is difficult to distinguish this scenario from the 
regular model discussed earlier because of the small dif- 
ferences in temporal/spectral behaviors of the GeV emis- 
sion. 



B. Other high energy processes in GRBs and 
afterglows 

1. High energy 'y—rays from pion production 

Katz [l60| | introduced pion production to explain the 
delayed 18 GeV photons of GRB 940217. In his model, 
these energetic photons resulted from collisions of rela- 
tivistic nucleons with T ^ 300 with a dense surrounding 
cloud, producing tt*', and decaying to tens GeV photons 
[l6^. 

The GRB outflow might be neutron-rich [161]. In a 
neutron-rich fireball, if the initial entropy > 400, the neu- 
trons and protons acquire a relative drift velocity causing 
inelastic n, p c ollis i ons a nd creating 7r° and yielding ~ 10 
GeV photons p^ . fT63j . 

In GRB internal shocks, protons can be accelerated to 



very high energies and photonpion collision create tt" or 
A"*" which in turn produces high energy photons. These 
photons are so energetic that they can not escape from 
the firebaU (see eq.dSO])) Dermer & Atoyan [Hi] 

argued that the ultra-high energy neutrons created in 
the p -|- 7 process play a crucial role in producing GeV 
emission. These neutrons are not confined by the mag- 
netic field of the blast-wave shell and flow out, and are 
subject to further photo-pion processes with photons in 
the surrounding medium to form charged and neutral pi- 
ous. The charged pious decay into ultra-relativistic elec- 
trons and neutrinos, whereas the decay of 7r° produces 
two 7— rays that are promptly converted into electron- 
positron pairs on the assumed Gauss-strength magnetic 
fields surrounding GRB sources. The synchrotron radi- 
ation of these energetic pairs can give rise to a strong 
GeV emission and may be able to account for the hard 
7-r ay c omponent detected in GRB 941017 jl64{ . loka et 
al. jl65j ] considered the /3-decay of the ultra-relativistic 
neutral beam (n —^ p + + D,,) and predicted extended 
GeV-TeV emission surrounding Gamma-Ray Burst rem- 
nants. 



2. Electromagnetic cascade of Te V 7— rays 

The electromagnetic cascade of TeV 7— rays was intro- 
duced by Plaga [73:3 as a model for the long-lasting MeV - 
GeV afterglow emission of GRB 940217 (see also [HI)- 
In this model TeV emission from the bursts c asca de on 
the CMBR. Dai & Lu (zj and Wang et al. [H3 pro- 
posed tha t the SSC emission of internal shocks (see also 
[l28l Il68j ]') or the very early forward shock may peak at 
TeV energies and then calculated the delayed MeV-GeV 
emission resulting in the electromagnetic cascade. De- 
tailed modeling of the electromagnetic cascade of GRB 
TeV emission has been carried out in [t^, [zl] . However 
this model requires a very low IGM field < 10""'* G, oth- 
erwise the resulting emission will be extremely weak (see 
eq.lpl). 



V. INTERPRETATIONS OF CURRENT 
OBSERVATIONS 

So far there are only few high energy GRB observations 
and in most cases the data is insufficient to carry out a 
detailed analysis. Most cases are consistent with SSC 
or EIC interpretation either from the prompt emission 
or from the afterglow shocks or from both. However, 
numerous other models have been put forwards. 

• Gonzalez et al. [2^ discovered a significant sub- 
GeV emission in 26 bright GRBs (including GRB 
941017) during the prompt phase. In these cases 
the high energy spectra are consistent with the 
single power law component observed by BATSE. 
The simplest and the most natural interpretation 
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FIG. 12: The SSC radiation of the forward shock undergoing 
energy injection and with evolving shock parameters, the case 
of GRB 940217: the thick solid hne is the light curve and the 
inserted plot is the spectrum (the times have been marked 
in th e plot), in the energy range of 30 MeV - 30 GeV (from 





FIG. 13: The SSC and the second order inverse Compton 
emission of internal shocks that have a synchrotron radiation 
peaking in UV/optical band (from [H^). Note that in this 
figure the photons' energy, in units of rrieC? ^ is measured in 
the rest frame of the GRB ejecta and not in the observer's 
frame. 



of the prompt high energy 7— rays is an internal 
shock synchrotron radiation in the MeV-GeV en- 
ergy range. Similar conclusion can be drawn for 
the GeV excesses seen in GRB 920622 and GRB 

940301 ;il. 

• The discovery of a 18 GeV photon about one hour 
after GRB 940217 is very amazing. The other 17 
high energy photons detected in phase-2 have sim- 
ilar energy ~ 100 MeV (see FigI3]). Excluding the 
18 GeV photon, the count rate and the energy of 
the high energy photons are almost constant. 

A slowly decaying MeV— GeV afterglow light curve 
is possible in the standard afterglow shock model 
if Vc.ssc < V < I'm, SSC (sec Table |V| . This inter- 
pretation, however, suggests (3 ~ 1/2, which is in- 
consistent with the o bserv ation /3 ~ 1.8 and thus is 
unlikely. Wei & Fan [l47l | suggest that the SSC of 
a modified forward shock can reproduce the data. 
Their numerical results are shown in Fig [121 An 
alternative is that the MeV-GeV plateau was the 
SSC radiation component of an X-ray plateau like 
that detected in GRB 070110. The two models dif- 
fer in the origin of the seed photons. In one case 
it is the forward shock X-ray emission and in the 
other a central engine afterglow. 

• In GRB 941017, both the duration and the fluence 
of the high energy emission component are about 
3 times that of prompt soft 7— rays. The spectrum 
of these hard gamma-rays is unusually hard {F^^ cx 

I/O). 

Granot & Guetta [l35| | (see also [l37| ) suggested 



that inverse Compton from reverse shock is most 
likely to provide the right temporal behavior. Be- 
cause of the very large energy emitted in hard 
7— rays, the reverse shock has to be relativistic. 
The second possibility is the inverse Compton from 
very early forward shock while the reverse shock 
is still going on. In both cases, the prompt soft 
7— rays overlap the shocked regions (as seen in the 
data, it is indeed the case), and plays an impor- 
tant role by providing the seed photons for an EIC 
pro cess t hat gives rise to the strong GeV emission 
16^ . Il24| . Two additional advantages of the EIC 
model are: (a) As shown in eq. (pD)) and in FigfTTl 
the duration of the high energy emission is natu- 
rally much longer than that of the seed photons 
(here the prompt soft 7— ray emission); (b) The 
spectrum could be as hard as Fi, cx :/-i/2~i/3^ as 
shown in figure 13 of [2^ . 



Dermer & Atoyan 164l | proposed a neutral beam 
model (see section HVB 1[) to explain the very hard 
spectrum Fi, (x , in which a Gauss-strength mag- 
netic fields of the circumburst medium is needed. 

• Although no very high energy emission was directly 
detected from GRB 080319B the unique prompt 
spectrum of this bursts suggests that it has been 
accompanied by a very strong GeV emission that 
would have carried at least 10 times more energy 
in the GeV component than in the prompt MeV 
emission (see pL69j for details). 

GRB 080319B was located at redshift z = 0.937. 
The peak energy of the spectrum was E-p ^ 600 
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keV, and the photon indexes below an d ab ove Ep 

are 0.8 and 3.4 respectively ^7^. The 

peak prompt V-band {v^ ^ 5 x 10^^ Hz) emission 
reached ~ 5th magnitude [1705. The spectrum of 
the optical to soft 7— ray emission cannot be inter- 
preted as a simple synchrotron spectrum as the low 
energy optical component is much too large. 

A possible solution is that the UV/optical emission 
is the synchrotron radiation of the internal shocks 
while the soft 7— rays are the corresponding SSC 
emission, i.e, hvm ^ 10 eV and ft.fm.ssc ^ E-p ^ 600 
keV. So the typical Lorentz factor of the emitting 
electrons is 7e „, ~ {ym,ssc/vmY/'^ ^ 200. The 
unusual spectrum of the 7— rays suggests a cool- 
ing Lorentz factor 7^^^ '-^ 7e,m- On the other 
hand, 7;^ ~ 500 T%^R^QL:^\^^^[e{\ + r2„dic)']-^ 
To derive this relation, the term (1 -I- Vssc) of 
eg. dTZ)) should be replaced by (1 -I- F2ndic) be- 
cause in this particular burst the 2nd-order IC is 
very important, as shown below. Given the ob- 
servation hv^F^^ / EpFs^ ~ 1/100 and given that 
Fjj cx v^/'^ at energies < /ii^m, we have Y^sc ^ 
100/(t'm/i^v)^^^ ~ 10- Note that these estimates 
are not very sensitive to the choice of the (un- 
known) value of Vm . Correspondingly e — s-b/ £c ^ 

■^2iidic ^&sc ^ 0.001 (which suggests, inciden- 
tally, that the outflow is not Poynting-flux domi- 
nated). A 7g J. ~ 200 thus requires that (Fq, R) ^ 
(10^, 10^'^ cm). The 2nd-order IC scattering of the 
internal shock electrons is still in the Thompson 
regime since ^'^^Ep/To < rricC^. The spectrum 

peaks at ~ 27'g ^Ep ~ 40 GeV and can escape the 
fireball freely. The luminosity of the 2nd order IC 
emission may be as high as ~ 10^^ erg s~^. Accord- 
ing to this model a significant GeV emission might 
have been detected by AGI LE i f it was not occulted 
by Earth at that moment jl69l | . It is very interest- 
ing to note that such a scenario has been outlined 
by Stern & Poutanen [6§|, as shown in FigfT3l 

Two additional high energy components with a du- 
ration ^ 100 s (i.e., longer than duration of the 
prompt 2nd-order IC emission) could arise in this 
burst. The optical emission in the time range of 
100 — 1000 sec is probably the high latitude emis- 
sion of the reverse shock (170i] . This interpretation 
implies that the reverse shock was relativistic and 
there was an overlap between the prompt emission 
and the reverse shock region. In this case, the re- 
verse shock electrons would be effectively cooled by 
the prompt photons (see section lTV A 2[) and would 
give rise to GeV-TeV EIC emission with a luminos- 
ity ~ 10^^ erg s~^. The huge amount of prompt 
UV/optical photons would cool the forward shock 
electrons effectively and give rise to another GeV- 
TeV EIC emission with a luminosity ~ 10^^ erg s~^ 
(see [HI). 



VI. SUMMARY AND OUTLOOK 

We have reviewed the observation of the prompt and 
afterglow high energy emission from GRBs (Sec. |IT]), then 
concentrated on the possible physical processes giving 
rise to these signals (Sec. IIII[) and applied them to GRBs 
and their afterglows (Sec. IIV[) . The likely interpretations 
of the observations have been presented in Sec. |Vl 

After reviewing various possible sources for very high 
energy emission we find, somewhat expectedly, that the 
most important source for very high energy emission in 
GRBs and their afterglows is inverse Compton either in 
the form of SSC or EIC. Our discussion has been based 
on the standard baryon-rich fireball model, which works 
pretty well in modeling late time GRB afterglow data 
(see d, [nl, [l3] for reviews) . 

The high energy emission is expected to be strong in 
the prompt emission phase and the early afterglow. The 
current limited data can be interpreted within the frame- 
work of the fireball model, though at times some modifi- 
cations are needed. Future high energy observation, par- 
ticularly in the early afterglow phase, will impose tighter 
constraint on these modifications. 

We expect that the internal shock SSC emission 
in both prompt emission phase and in central engine 
afterglows will be detectable by LAT in the MeV-GeV 
energy range (see Table IIIip . The higher than 50GeV 
SSC emission from internal shocks is expected to be 
detectable only for some extremely bright GRBs with a 
Fo > 10^ and a z < 1. A rough estimate of the detectable 
count rates of high energy photons from GRB forward 
shock by LAT and MAG IC is summarized in Table |Vll 
Recently Xue et al. [l7l| calculated the VHE SSC emis- 
sion of the forward shock of specific nearby GRBs having 
a z < 0.25 and found that a significant detection was 
expected only in GRB 030329. The results are consistent 
with the null detection of MAGIC [H, Whipple [U and 
H.E.S.S. m. MAGIC-II and H.E.S.S.-II wiU lower the 
energy threshold to about 30 GeV. With a very large 
effective area ~ (1 — 5) x 10^ cm^ and the much less 
absorption by the IR-background, significant detections 
of the tens GeV photons from GRBs and afterglows will 
be possible. We may be able to use these detections to 
calibrate the IBL models at high redshifts. Below we 
focus on the MeV-GeV signatures that are detectable by 
LAT and the possible constraints on the astrophysical 
model that might arise from this new data. We divide 
this discussion to the prompt emission and the afterglow: 

Prompt Emission 

The field of view of GBM is all sky not occulted by the 
earth and that of LAT is ~ 2.5 sr. So - 1/5 GRBs wifl 
be within the field of view of both GBM and LAT. With 
a good quality keV— GeV spectrum, we may achieve the 
following goals: 

• The particle acceleration process in GRBs can be 
better constrained. In principle, the electrons can 
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TABLE VI: Expected signals from GRB forward shock SSC emission. The absorption of the VHE photons by infrared 
background has been taken into account, based on Table|IIl These values are calculated for a typical burst with i5k ~ 10^^ erg (at 
the end of the X-ray shallow decline) and z = 1, correspondingly to the burst with a sub-MeV 7— ray fluence of ~ 10~^ erg cm~^ 
(from 24]). Note that we have used for Magic an effective area, in the energy range of 50GeV-10TeV, of Sdot ~ 10^ m^ This 
might be an overestimate for ~ 50 GeV. 



iVdct(> 20MeV) iVdct(> lOGeV) iVdct(> 50GeV) 
LAT LAT (Im^) MAGIC (lO^m^) 



Standard afterglow: ISM 


~ 13 


~ 10 X 10-^ 


~ 11 


Standard afterglow: wind 


~ 12 


~ 8 X 10^2 


~ 7 


Energy injection: ISM 


~ 3 


~ 2 X 10-2 


~ 2 


Energy injection: wind 


4 


~ 2 X 10-2 


- 2 


Time increasing Sc'- ISM 


~ 3 


~ 2 X 10-2 


~ 2 


Time increasing £c: wind 


~ 3 


~ 1 X 10^2 


~ 1 



be accelerated to an energy ^ 

205/-v^ TcV (mea- 
sured in the rest frame of the emitting region) and 
can give rise to GeV synchrotron radiation (see 
eq.(16)). A cut-off at much lower energy (say, a few 
MeV) will challenge the internal shock model or the 
currently adopted particle acceleration model, that 
assumes dn' /d^'^ cx 7e~^. 

• The escape of the most energetic prompt photons 
from the fireball depends strongly on To the bulk 
Lorentz factor of the GRB ejecta (see eg. ljSU)) ) j82j 
and on t he radius of the internal energy dissipa- 
tion [l72j |. A detected GeV cutoff would constrain 
these two parameters (unless the redshift is suffi- 
ciently large that the cosmological attenuation will 
give rise to such a cutoff too (see Figl^])). 

• A strong SSC GeV emission component is expected 
in the standard baryonic dominated model. A lack 
of such emission in the prompt spectrum may favor 
the internal magnetic energy dissipation model for 
the prompt 7— rays [11,111, in which the inverse 
Gompton effect is weak (unless Fq or dt^ is very 
small (see eg. ([50)) ). 

Actually, as shown in the last paragraph of section 
El the prompt GeV emission of GRB 0803 19B may 
have a luminosity as large as ~ 10^'^ erg s^"'^ or even 
larger, which is '--^ 10 times more energetic than the 
prompt optical— MeV emission. AGILE GRID should 
have detected such an amazing component if it was not 
occulted by the earth at that moment. Such a detection 
would be crucial to pin down the prompt emission 
mechanism (s). 

Afterglow 

The afterglow, and in particular the early afterglow is 
an additional powerful source of very high energy emis- 
sion. Motivated by the identification of a cano nical X-ray 
afterglow light curve of Swift GRBs [loi,[io3. Fan et al. 
[2^ argued that a similar high energy emission light curve 
should be observed (FigfTi)). Given the small number of 
expected high energy photons (typically ~ 10 for a bright 
burst at redshift z ^ 1) detectable for LAT, these novel 



features are not likely to be identified as frequently as in 
X-ray band. However, a detection of more than 10^ — 10"^ 
sub— GeV photons is possible in some extremely bright 
bursts and these case can be used to test the prediction 
shown in FiefHl 

• EGRET had discovered long-lasting MeV-GeV 
plateau in GRB 940217 and an energetic delayed 
sub-GeV plateau in GRB 941017. These are good 
candidates of a shallow decline (or plateau) phase 
and very early EIC plateau of the high energy after- 
glow light curve shown in Fig[T4l LAT is expected 
to detect more such cases in the coming years. 

• The SSG emission of central engine afterglows, such 
as X-ray flares and X-ray plateaus followed by a 
sharp drop, may give rise to detectable sub-GeV 
flares and plateaus. The EIG emission component 
usually lasts much longer than the seed photon 
pulse and may be outshined by the SSC emission 
component of the forward shock (see Fig lTT]) . As 
long as Yssc > 1 and the forward shock electrons are 
in fast cooling, the EIC process can not enhance the 
detectability of the high energy afterglow a lot be- 
cause in the absence of EIC, the SSC will radiate 
significant part of the energy Lein into GeV ener- 
gies. Distinguished EIC signatures are expected if 
the forward shock cools very inefficiently before the 
flare. 

• The establishment of a canonical high energy af- 
terglow light curve in some extremely bright GRBs 
will confirm current interpretations of the peculiar 
Swift X-ray afterglow data. 

The prospects for these advanced are good. The very 
high emission might provide essential clues to the na- 
ture of GRBs and in particular to the elusive conditions 
within the emitting regions and the emission processes. 
However, as suggested by past experience new challenges 
and surprises are bound to emerge when a new observa- 
tional window is opened and new observations become 
available. GLAST LAT, on one hand and new improved 
Cenrenkov telescopes, with a larger collection area and 
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FIG. 14: The expected high energy emission sig natures are shown in the lower panel, corresponding to the schematic X-ray 
light curve shown in the upper panel (from l24|). Please note that the SSC flares could be in GeV-TeV energy range if the 
X-ray flares are powered by external shocks. 



a lower energy detection threshold, on the other, may 
provide such surprises in the near future. 
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